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i.  Drrtorocnow 


overview  or  nmnaoiTAL  approach  and  aqoiwn 

— Four  ^  distinct''  experiments  art  performed  to  determine  the  mechanical 
properties  of  lur.g  tissue: 

1*  Measurement  of  the  incremental  bulk  Modulus ,  K.  The  equlpisent  Imposes 
Incremental  changes  to  the  lung  volume  and  measures  the  corresponding 
transpulaonary  pressure.  From  the  pressure-volume  curves  the  bulk 
modulus  Is  computed.^ Figure  1-1  shows  a  photograph  of  this  equipment. 
Details  of  analysis  and  results  are  presented  in  Section  2.1. 

2.  Measurement  of  the  incremental  Young's  modulus,  E.  Two  types  of  mea¬ 
surements  are  done: 

a.  A  known  vertical  deflection  is  applied  on  the  pleural  surface  and 
the  corresponding  vortical  load  Is  measured. 

b.  A  circular  flexible  polyethylene  membrane  of  various  diameters 
(7  mm  to  16  am,  depending  on  the  also  of  the  lung)  Is  cemented  to 
a  flat  portion  of  the  pleura,  and  a  vertical  deflection  Is 
applied  on  the  center  of  the  membrane^)  The  corresponding  vertical 

— load  Is  measured.  Figure  1-^2  shows  a  photograph  of  this  equip¬ 
ment.  Performance  has  been  thoroughly  tested  and  Is  satisfactory. 
Details  of  the  experiments  and  their  results  are  presented  In 
Section  2.2. 

3.  Measurement  of  material  constants  for  human  lungs  in  large  deforma¬ 
tion.  This  consists  of  two  parts: 

a.  The  TRIAX  testing  machine,  designed  and  constructed  in  our  labor¬ 
atory  earlier,  was  used  to  measure  the  tissue  components  of  the 
human  lung  stress-strain  relationship i,  FI guv a  1-3  shows  the 
equipment;' "TWW 'effort itas  made  to  adapt  a  new  computer  In  the 
laboratory,  PDP  1103,  to  perform  on-line  experiments  and  analy¬ 
sis.  ThlB  has  been  accomplished.  Results  are  contained  In  Section 
2.3. 

— i  b.  Pressure-volume  curves  of  air-filled  human  lungs  are  measured. 

This  work  is  done  In  Dr.  Scbin's  laboratory  In  Los  Angeles  where 
fresh  specimens  can  be  obtained  and  tested  as  soon  as  they  arrive 
In  the  laboratory.  ,  , 
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4.  Measurement  of  Wavs  Propagation  In  tha  Lung .  For  convenience  of  de- 
tailad  examination,  axclaad  lunga  ara  used.  k  small  shock  tuba  ia 
uaad  to  laposa  Impact  loading  on  tha  lung.  Prasaura  transducara  ara 
uaad  to  Sanaa  tha  ahock  and  tha  atraaa  wavas.  Tha  spaad  of  wavs  prop** 
agatlon  la  dataralnad  from  tha  first  arrival  tlaa  of  tha  prassura 
wava  front  and  known  diatanca  batwaan  saneora.  Figure  1-4  shows  ona 
of  tha  two  shock  tubas  daaignad  and  constructed.  Flgura  1-5  ahows  tha 
tast  setup  for  tha  lung.  Results  ara  presented  In  Section  3, 


Figure  1-2.  Kqulpaant  for  Maturing  IncreMntal  Toung't  nodulut  of  lung 
parenchjM  by  lndantatlon  and  by  flaxlbla  dlak 
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2.  UJRG  PARENCHYMA  MECHANICAL  PROPERTIES 


2.1  INCREMENTAL  BULK  MODULI  OF  ELASTICITY  OF  THE  LONG 
2.1.1.  Definitions 

In  the  lung,  the  stress  acting  on  any  cross  section  consists  of  two 
parts:  the  stress  due  to  stretching  of  the  lung  tissue  (°ti8sue)>  and  the 
alveolar  gas  pressure  (F^).  We  write 

stress  -  otifl8ue  -  a?A  ,  (2.1-1) 


where  a  Is  the  fraction  of  the  cross  section  that  Is  exposed  to  the  alveolar 
gas,  1-a  Is  the  fraction  occupied  by  the  Interalveolar  septa.  See  Figure 
2.1-1.  On  the  boundary  of  the  lung  is  the  visceral  pleura.  One  side  of  the 
pleura  Is  exposed  to  the  pleural  pressure  pPL.  The  other  side  Is  exposed  to 
the  stress  given  by  Eq.  (2.1-1).  The  tension  in  the  pleura  can  be  expressed  In 
terms  of  the  principal  stress  resultants,  Tj  and  T2  per  unit  length.  Let  the 
principal  curvatures  of  the  pleura  be  1/Rj,  and  1/R2,  respectively.  Then  the 
condition  of  equilibrium  of  the  pleura  Is  (Fig.  2.1-2): 


°ti8sue  apA  +  PPL 


(2.1-2) 


Due  to  the  weight  of  the  lung,  or  motion  of  the  chest,  the  stress  in  the  lung 
Is  nonuniform  and  the  equilibrium  of  the  lung  requires  the  tension  in  the 
pleura.  But  If  the  right-hand  side  of  Eq.  (2.1-2)  vanishes,  then 


°tl88ue  *  apA  PPL 


(2.1-3) 


°tlssue  *8  a  *unct:ion  °f  the  strains  in  the  lung.  If  homogeneous  strain  Is 
assumed ,  then  ot^gsue  la  a  function  of  the  lung  volume;  so  are  pA  and  pp^.  The 
volumetric  strain  of  the  lung  being  the  ratio  of  change  of  volume  of  the  lung 
to  the  lung  volume,  Aviung^Vlung ’  t^e  incremental  volumetric  (or  bulk)  modulus 
of  elasticity  of  the  lung  tissue  is  defined  as 


»  m  v  ^tissue 

tissue  lung  d^ung 


d(“PA  “  PPL) 


lung 


dV 


lung 


(2.1-4a) 


Figure  2.1-1.  Concept  of  macroscopic  stress  defined  over  a 
small  cross  section  of  the  lung  tissue 
enclosing  a  number  of  alveoli  and  ducts. 


Figure  2.1-2.  An  enlargement  of  a  small  element  of  pleura,  showing  exposure  to 
pleural  pressure  (ppL)  on  one  side,  and  to  alveolar  gas  (p^)  and 
Interalveolar  septa  on  the  other.  The  tension  In  the  pleura 
multiplied  by  the  curvature  of  the  pleura  Is  equivalent  to  a 
lateral  load,  and  contributes  to  the  balance  of  forces  In  the 
direction  normal  to  the  pleural  surface,  resulting  In  the  terms 
on  the  right-hand  side  of  Eq.  (2.1-2). 
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This  is  the  modulus  relevant  to  the  study  of  breathing,  and  is  normally  mea¬ 
sured  by  experiments  in  which  either  pA  or  ppL  is  varied,  with  the  airway 
open. 


■"n  the  study  of  response  of  the  lung  to  shock  waves,  it  is  relevant  tc 
consider  the  case  in  which  the  glottis  is  closed.  The  shock  process  is  so  fast 
that  the  bulk  movement  of  gas  through  the  glottis  is  negligible  in  the  first 
few  moments.  In  the  closed  glottis  experiments,  we  may  define  the  following 
Incremental  volumetric  moduli  of  elasticity: 


K. 

lung 


dP, 


lung  dV 


lung 


(2.1-4b) 


*PL- 


dp 


PL 


lung  dV. 


lung 


(2. l-4c) 


We  shall  call  Kp^  the  "chest  bulk  modulus."  KpL  is  of  interest  to  the 
calculation  of  the  chest  volume  in  response  to  shock  waves  enveloping  the 
chest.  Ktiasue  la  0f  interest  to  the  assessment  of  the  stresses  in  the  lung 
tissue,  and  hence  of  Its  strength  and  damage. 


From  Eqs.  (2.1-4a),  (2.1-4b),  and  (2.1-ac)  we  obtain 

aK  +  K__  -  K  ,  . 

lung  PL  tissue 


(2.1-5) 


2.1.2.  Open  Airway  Experiments 

An  experiment  was  designed  to  impose  increments  of  volume  to  the  lung 
while  allowing  air  transfer  between  the  inside  and  outside  of  the  lung.  An 
excised  lung  is  suspended  inside  a  lucite  box.  Air  volume  in  the  lung  is 
varied  by  a  pump  connected  to  the  trachea  cannula  while  airway  pressure  (pA) 
is  measured  and  external  pressure  (Pp^)  Is  held  atmospheric.  The  bulk  modulus 
is  computed  from  pressure-volume  curves  obtained  in  this  fashion.  Figure  2.1-3 
shows  a  typical  p-v  curve. 

Results  are  shown  in  Table  2.1-1.  It  is  seen  that  the  magnitude  of  the 
bulk  modulus  of  the  lung  tissue  of  the  rabbit  is  of  the  order  of  3  or  4  times 
the  transpulmonary  pressure.  Lai-Fook  (1979)  has  shown  that  Ktissue  dog 
lung  Is  of  a  similar  magnitude.  Table  2.1-2  shows  the  results  for  cats  and 
humans . 
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AP  TRANSPULMONARY  PRESSURE 
cm  HjO 


Figure  2* 1-3.  Frasaure-voluae  curvet  of  air-filled  human  lung 


Figure  2.1-4.  Schematic  drawing  of  the  hulk  modulus  measuring  equipment. 


Table  2.1-1*  The  Bulk  Modulus  of  the  hung  Tissue  (Ktisgus)  of  the  Rabbit 
Determined  by  Open  Airway  Experiments 


Transpulaonary 

Pressure 

PA  "  PPL 
(cm  H2O) 

n 

^tissue  <cm  H2°> 

Mean 

±  S.D. 

Inflation 

4 

4 

13.60 

2.42 

Deflation 

4 

4 

12.26 

0.65 

Inflation 

6 

4 

21.32 

1.62 

Deflation 

6 

4 

15.07 

0.46 

Inflation 

8 

4 

30.69 

2.21 

Deflation 

8 

4 

28.35 

1.31 

Inflation 

10 

4 

36.87 

3.70 

Deflation 

10 

4 

40.99 

2.15 

12 


Table  2.1-2.  The  Balk  Modulue  of  Lang  Tissue 


e. 

Cats'  Lung  Tissue 

franspulao«.*ry 

Pressure 

^tissue 

(cm  R2O) 

_ FA  -  fPL  <“  V> 

n  Dean 

S.  D. 

Inflation 

4 

3 

♦ 

15.62 

2.45 

Deflation 

4 

3 

12.19 

2.07 

Inflation 

6 

3 

32.27 

3.83 

Deflation 

6 

3 

26.16 

5.15 

Inflation 

8 

3 

49.59 

7.38 

Deflation 

8 

3 

46.03 

11.52 

Inflation 

10 

3 

64.47 

12.13 

Deflation 

10 

3 

61.89 

15.30 

b.  Human 

Lung  Tissue 

Transpulmonary 

Pressure 

'A  -  PPL  H2°> 

HjO) 

Inflation 

5 

18.10 

Deflation 

5 

21.94 

Inflation 

6.5 

21.77 

Deflation 

6.5 

25.40 

Inflation 

8 

26.19 

Deflation 

8 

28.73 

Inflation 

9.5 

29.57 

Deflation 

9.5 

29.80 

I.  CLo— J  Airway  hp>rimti 
Iqnlp— at  amd  Method 

The  equipment  sketched  in  Figure  2.1-3  was  used.  The  volume  of  sir  in  the 
box  end  external  to  the  lung  can  be  varied  by  a  Harvard  Infusion-withdrawal 
pump. 

The  calculation  needs  *the  volume  of  the  lung  tissue.  This  was  obtained  by 
measuring  the  volume  of  a  completely  atelectatic  lung.  The  lungs  of  the  animal 
under  anesthesia  were  ventilated  with  100%  pure  oxygen  for  one  hour  and  then 
the  trachea  was  clamped  off.  The  gas  In  the  lungs  was  then  quickly  absorbed  by 
the  circulating  blood  and  the  lungs  collapsed  completely.  After  excision,  the 
volume  of  the  tissue  was  determined  by  Immersing  it  in  water  and  measuring  the 
rise  of  water  surface  in  a  tube. 

The  lung  was  then  hung  in  the  airtight  lucite  box.  A  known  volume  of  air 
was  then  Injected  using  a  syringe  so  as  to  open  all  airspaces,  creating  an 
airway  pressure  of  about  20-25  cm  ^0  while  the  pleural  pressure  was  atmo¬ 
spheric.  For  each  experiment,  a  known  volume  of  air  was  added  or  withdrawn 
from  the  lung  using  a  syringe  to  obtain  a  desired  initial  pA  while  ppL  »  0. 
The  airway  valve  was  then  closed.  Further  change  of  lung  volume  was  imposed  by 
the  infusion-withdrawal  pump.  In  each  change  of  volume  the  corresponding 
transpuluonary  pressure  (tracheal  pressure  minus  pleural  pressure)  and  the 
airway  pressure  were  me*  mreJ  simultaneously  using  the  Valldine  pressure 
transducers.  During  tlv.  experiment,  the  airway  pressure  was  recorded  using  a 
strip-chart  recorder,  and  the  pump  volume  was  plotted  against  pulmonary  pres¬ 
sure  using  an  x-y  jecorier.  An  x-y  recorder  was  used  also  to  plot  change  of  pA 
versus  charge  of  pa  -  pp^.  These  values  are  used  to  compute  the  bulk  moduli. 

Analysis 

The  volume  of  the  lung  la  the  sum  of  the  volume  of  gas  in  the  lung 
^alv.gas*  and  the  voluoiB  of  the  lun«  tissue  (vti88ue).  Hence 

V.  -  V  ,  +  V„. 

lung  alv.gas  tissue 

Since  the  airway  and  the  external  chamber  are  both  closed,  we  have,  by  gas 
law, 

pa  •  T.iv.g..  -  H*Rr  •  <2->-6> 

u 


Hence ,  when  pj^,  PA,  and  \u.gmm  •«»  changed  by  Ap^,  ApA,  end  ^valv>gM: 


have 


fp.  +  ApJfV  -  +  AV  ,  „  )  -  M.RT  ; 

rA'v  elv.gea  alv.ga* '  A 


we 


AV 


M  RT 

A 


-  V. 


lung  pA  +  ApA  elv.gae 


AP, 


PA  +  ApA  elv.gea  * 


(2.1-7) 


Substituting  this  Into  the  definition  of  cheat  bulk  modulus  Eq.  (2.1-Ac),  we 
obtain 


elv.gae  A 

On  the  other  hand,  for  lungs  with  closed  airways,  Eq.  (2.1-6)  laplles 


(2.1-8) 


dp.  dV  , 

— »  -  0  . 

PA  »lv.g«t 

Hence  on  substituting  Into  Eq.  (2.1-Ab),  we  obtain  the  bulk  modulus  of  the 
lung  In  the  case  of  closed  glottis: 

V. 


hmm  ’  -  T-f**-  **  • 

alv.gas 


Then  by  Eq,  (2.1-5), 


V,  *P„, 

m  lung  PL 


'tissue  V,_ _  Ap»  VKA  A 


(p.  +  ApJ  -  op.  . 


(2.1-10) 


alv.gas  A 

Since  o  Is  an  unknown  anatomical  parameter,  we  compute  a  K  for  closed  airway 
defined  by  Eq.  (2.1-10)  with  o  ■  1: 


K 


Vlung  A**FL  r  ,  r  1 

v~~ - JT  lp*  +  4V  '  pA 

alv.gas  A 


then 


K  -  K  +  (1  -  u)p.  - 

tissue  rA  V 


JLung. 


alv.gas 


(2.1-11) 


(2.1-12) 
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Table  2.1-3  show*  a  sat  of  typical  results  on  an  animal.  It  shows  the 
similarity,  with  opposite  signs,  of  the  values  of  the  chest  hulk  modulus  and 
lung  gas  modulus  when  the  airway  is  closed.  The  sum  yields  the  bulk  modulus  of 
the  lung  tissue. 

Table  2.1-4  summarises  the  experimental  reaults  of  the  rabbit  lung. 

Plscuasion 

The  chest  bulk  modulus,  RpL,  for  lung  experiments  with  glottis  closed,  is 
seen  to  be  much  larger  than  the  bulk  modulus  of  the  lung  tissue, 
is  nearly  equal  to  the  bulk  modulus  of  the  air  In  the  lung. 

The  values  obtained  from  closed  airway  experiments  are  seen  to  be 

different  from  those  obtained  from  open  airway  axperimanta.  The  reason  for 
this  discrepancy  Is  not  clear.  The  major  difference  to  the  lung  tissue  In 
these  two  experiments  Is  that  when  the  airway  la  closed  very  large  changes  In 
both  pleural  and  airway  preasuras  are  required  to  obtain  a  small  change  In  the 
transmural  pressure,  pA  -  Pp^.  whereas  in  the  open  airway  experiments  is 
kept  constant  and  ApA  is  exactly  equal  to  A(p^  -  ppL).  Lung  inflation  is,  as 
is  seen  from  Kqs,  (2.1-2)  and  (2.1-3),  essentially  a  function  of  the  trans- 
mural  pressure  PA  -  ppL,  but  blood  vessels  respond  to  (blood  pressure  -  pA)  in 
the  capillaries  and  to  (blood  pressure  -  pp^)  In  the  large  vessels.  Hence  the 
blood  vessels  wist  be  deformed  to  a  very  different  degree  In  these  two  cases. 
The  pleura  contains  bronchial  and  lymph  vessels  and  must  recognise  the  differ¬ 
ence  similarly.  Further,  we  have  lumped  the  airway  (the  part  In  the  test  box) 
with  the  lung  In  the  data  analysis ,  and  treated  the  stress  in  the  lung  as 
homogeneous.  Whether  these  factors  explain  the  difference  in  *ti88ue  from  open 
and  closed  airway  experiments  or  not  remains  to  be  investigated. 

2.2  DRUniAL  YOUWG'S  MODULUS  OF  ELASTICITY 

Once  the  bulk  modulus  K  of  the  lung  tissue  is  determined,  it  is  possible 
to  obtain  results  for  the  shear  modulus  G,  Young's  modulus  E,  and  Poisson's 
ratio  v  through  the  Indentation  test. 


2.2.1.  Apparatus 

A  photograph  of  the  equipment  used  in  this  test  is  shown  in  Figure  2.2-1. 
It  consists  of  a  Mitutoyo  deflection  apparatus  (accurate  to  a  0.005  tom  deflec¬ 
tion),  Statham  pressure  transducer  (Model  UC3),  a  Hewlett  Packard  31 1A  trans¬ 
ducer  amplifier-indicator,  a  Hewlett  Packard  7100  BM  strip  chart  recorder,  and 
a  manometer,  which  was  used  to  inflate  the  lung  to  various  pressures.  The 
pressure  transducer  had  a  highly  sensitive  (3.175  mm  diameter)  probe  which  was 
mechanically  attached  to  the  deflection  apparatus.  This  setup  made  it  possible 
to  measure  the  resistive  load  of  the  lung  tissue  resulting  from  an  applied 
(known)  deflection.  The  signal  from  the  pressure  transducer  was  then  amplified 
and  transmitted  to  the  strip  chart  recorder  where  the  load  could  be  accurately 
measured . 

The  pressure  transducer  was  calibrated  with  a  dual-arm  balance;  a  50  gm 
weight  was  placed  on  one  end  while  the  probe  was  placed  on  the  otl. .  .*  end.  The 
deflection  apparatus  was  utilized  to  move  the  probe  up  or  down  until  the  bal¬ 
ance  arm  was  completely  leveled,  which  implied  that  the  probe  was  sensing  the 
full  50  gm.  The  range  of  the  strip  chart  recorder  was  then  reduced  since  only 
small  deformations  were  involved.  Under  such  conditions  the  stress-strain 
relationship  of  the  lung  is  considered  to  be  linear. 

2.2.2.  Experimental  Method 

In  each  experiment.,  a  fresh  lung  was  obtained  with  the  airway  completely 
opened.  The  manometer  was  used  to  inflate  the  lung  to  the  desired  pressure. 
The  lung  was  placed  firmly  on  a  tray  so  that  when  a  deflection  was  applied  to 
the  lung,  the  corresponding  measured  resistive  load  would  be  due  entirely  to 
the  tissue.  The  lung  was  kept  continuously  wet  with  saline  to  preserve  its 
freshness.  However,  the  saline  was  removed  from  the  tray  prior  to  applying  a 
deflection  (to  the  tissue)  since  water  trapped  beneath  the  lung  would  result 
in  an  incorrect  load  measurement.  A  polyethylene  disk  of  (known)  diameter  d 
and  negligible  thickness  was  glued  to  the  flat  surface  of  the  lung  using  Kodak 
Eastman  910  adhesive.  Though  only  small  deflection  was  applied  to  the  lung, 
the  disk  was  used  to  protect  the  lung  tissue  from  being  damaged  by  the  probe, 
which  would  also  result  in  faulty  load  measurements.  The  disk  and  adhesive 
used  did  not  appreciably  affect  the  elasticity  of  the  lung  tissue. 


The  lung  was  placed  underneath  the  deflection  apparatus  with  the  probe 
situated  directly  over  the  center  of  the  disk.  The  reference  point  was  deter¬ 
mined  by  applying  a  downward  deflection  until  the  transducer  just  sensed  a 
positive  load.  The  deflection  apparatus  was  utilized  to  move  the  probe  back  up 
to  the  point  where  the  recorder  indicated  zero  load.  This  was  the  probe's  ref¬ 
erence  point:  any  positive  deflection  would  result  in  a  positive  load.  A  known 
deflection  (w  in  mm)  was  then  applied  to  the  center  of  the  disk  and  the  corre¬ 
sponding  load  (L  in  gm)  was  obtained  from  the  recorder  after  a  static  equilib¬ 
rium  condition  was  reached.  To  obtain  consistent  results,  the  lung  was  precon¬ 
ditioned  before  applying  each  deflection.  This  also  eliminated  any  collapsed 
lung  tissue  beneath  the  disk  that  was  created  by  previous  deflections.  (As 
mentioned,  the  disk  greatly  reduced  this  effect.)  For  each  lung,  tests  were 
run  at  several  constant  Inflation  pressures  while  the  pleural  pressure  was 
kept  atmospheric. 

2.2.3.  Data  Analysis 

Knowing  the  bulk  modulus  K,  the  other  incremental  moduli  were  then  calcu¬ 
lated  uoing  the  following  equations 

v  -  A  -  (A2  -  A  +  *)1/2  A  -  3  Kwdit/IOOL  (2.2-1) 

E  -  6K(0.5  -  v)  (2.2-2) 

G  -  E/2(l  +  v)  (2.2-3) 

The  indentation  test  was  performed  on  several  rabbit  and  cat  lungs,  and  on  a 
single  human  lung.  Typical  plots  of  the  indentation  test  of  L  versus  w  are 
shown  in  Figures  2.2-2  and  2.2-3.  They  are  all  extremely  linear,  as  theory 
predicts  for  infinitesimal  deformations  (Hooke's  Law).  The  values  of  w/L  in 
Eq.  (2.2-1)  were  obtained  by  the  inverse  slopes  of  the  least  squares  equation 
given  next  to  each  plot.  As  expected,  the  slopes  of  these  plots  become  steeper 
at  higher  inflation  pressures  because  the  lung  becomes  stiffer,  and  for  a 
given  deflection,  the  resulting  resistive  load  of  the  tissue  will  be  larger. 

An  average  of  L/w  was  obtained  at  each  inflation  pressure  from  all  of  the 
tests  performed.  This  mean  value  was  then  used  in  Eq.  (2.2-1)  for  the  calcula¬ 
tion  of  v  and  subsequently  for  E  and  G. 
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6  cmHjO 


L- 1.222  W 


Pa-Ppl-  4  cm  H20 
L  *=  0.929  W 


0  0.5  1.0 

DEFLECTION,  W  (mm) 


Figure  2.2-2. 


Experimental  results  on  the  load  (L)  vs  deflection  (w)  obtained 
by  the  anal 1-indent at Ion  method  described  in  the  text.  Rabbit 
lungs.  Equations  obtained  by  least -squares  fitting  for  various 
transpulmonary  pressures  (pA  -  ppL)  are  listed  in  the  figure 
for  the  four  cases. 


LOAD  (9m) 


2.2.4.  Results 


Tables  2.2-1  through  2.2-3  list  the  calculated  mean  results  of  the  incre¬ 
mental  moduli  for  the  rabbit,  cat,  and  human  lung,  respectively.  Notice  that 
in  each  table  the  values  of  E  and  G  are  very  constant  for  each  Inflation  pres¬ 
sure  regardless  of  whether  the  inflationary  or  deflationary  value  of  K  is 
used.  Also,  the  values  of  E  and  G  both  Increase  with  the  inflation  pressure, 
which  is  expected  since  E  -  stress/strain  and,  as  previously  discussed,  for  a 
given  strain,  the  stress  will  be  larger  at  higher  inflation  pressures.  A  com¬ 
parison  of  the  tables  shows  that  the  values  of  Young's  modulus  at  each  infla¬ 
tion  pressure  are  similar  for  both  rabbit,  cat  and  goat  lungs.  Poisson's  ratio 
for  the  rabbit  lung  varies  from  0.318  to  0.400  and  has  a  general  tendency  of 
increasing  with  the  Inflation  pressure.  See  Figure  2.2-4.  The  values  of  v  for 
the  cat  lie  in  a  similar  range  (between  0.307  and  0.405)  and  have  a  definite 
pattern  of  increasing  with  the  inflation  pressure.  The  values  of  v  for  the 
human  lung  lie  in  a  smaller,  lower  range  between  0.276  and  0.333,  and  have  no 
such  definite  pattern.  In  fact,  those  values  of  Poisson's  ratio  for  the  human 
lung  corresponding  to  the  deflationary  values  of  the  bulk  modulus  actually 
decrease  with  increasing  inflation  pressure. 

2.3  LARGE  DEFORMATION  STRESS-STRAIN  RELATIONS 

Experiments  were  done  to  determine  the  mechanical  properties  of  human 
lung  parenchyma  under  large  deformation.  The  tissue  is  subjected  to  large 
deformation.  Rectangular  slabs  of  human  lung  tissue  3  x  3  x  0.4  cm  in  size 
were  cut  from  frozen  saline-filled  post  mortem  lungs  from  persons  17-41  years 
of  age.  The  slabs  were  subjected  to  biaxial  loading,  with  displacement  in  one 
direction  varying  sinusoidally  with  time  while  that  in  the  perpendicular 
direction  being  kept  constant.  The  resulting  deformations  were  monitored  con¬ 
tinuously  in  the  center  of  the  specimen  by  video  dimension  analyzers.  All 
force  and  displacement  data  were  recorded  digitally  on-line  by  a  PDP  11/03 
computer.  The  specimens  were  soaked  in  normal  saline  at  37°C  and  pH  6.3,  7.25, 
and  7.4,  The  resulting  stress-strain  curves  were  highly  non-linear.  Visco¬ 
elastic  features  of  hysteresis,  creep  and  relaxation  were  observed.  Hysteresis 
of  the  tissue  in  cyclic  loading  changed  only  to  a  minor  extent  when  the  strain 
rate  was  varied.  The  concept  of  quasllinear  viscoelasticity  and  pseudo- 
elasticity  (Fung,  1972)  seems  applicable  to  the  experimental  data  obtained.  An 


Table  2.2-1.  The  Incremental  Young's  Modulus,  Poisson's  Ratio  and 
Shear  Modulus  of  Rabbit  Lung  Parenchyma 


Trans- 

pulmonary 

Pressure 

PA  “  PPL 
(cm  H20) 

n 

E 

(cm  HjO) 

Poisson* 8 

Ratio 

V 

G 

(cm 

t 

\o) 

Mean 

S.D. 

Mean 

S.D. 

Mean 

S.D. 

Inflation 

A 

4 

12.719 

0.770 

0.363 

0.009 

4.666 

0.282 

Deflation 

4 

4 

12.794 

0.439 

0.338 

0.006 

4.781 

0.164 

Inflation 

6 

4 

16.659 

0.398 

0.376 

0.003 

6.053 

0.145 

Deflation 

6 

4 

16.883 

0.213 

0.318 

0.002 

6.405 

0.081 

Inflation 

8 

4 

20.527 

1.153 

0.376 

0.006 

7.459 

0.419 

Deflation 

8 

4 

20.551 

0.696 

0.371 

0.004 

7.495 

0.254 

Inflation 

10 

4 

22.882 

0.837 

0.389 

0.004 

8.237 

0.301 

Deflation 

10 

4 

22.819 

0.888 

0.400 

0.004 

8.150 

0.317 
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Table  2.2-2.  Incremental  Elastic  Moduli  for  Lung  Tissue  of  Cat 


Transpulmonary 


PA 

Pressure 
-  PpL  (cm  H2O) 

n 

E  (cm 
Mean 

h2o) 

S.D. 

v  (cm  H20) 
Mean  S.D. 

G  (cm  H20) 
Mean  S.D. 

Inflation 

4 

3 

13.71 

0.08 

0.351 

0.027 

9.26 

0.13 

Deflation 

4 

3 

13.84 

0.10 

0.307 

0.038 

9.04 

0.20 

Inflation 

6 

3 

19.00 

0.06 

0.401 

0.012 

13.31 

0.07 

Deflation 

6 

3 

19.13 

0.11 

0.375 

0.023 

13.15 

0.14 

Inflation 

8 

3 

27.78 

0.10 

0.405 

0.014 

19.52 

0.12 

Deflation 

8 

3 

27.85 

0.18 

0.395 

0.024 

19.43 

0.21 
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Table  2.2-3.  Incremental  Elaetlc  Moduli  for  Lung  Tlsaue  of  Human 


Trans pulmonary 
Pressure 

PA  -  PPI.  <=»  H2°> 

E 

(cm  H20) 

V 

(cm  H20) 

G 

(cm  H20) 

Inflation 

5 

22.15 

0.296 

8.55 

Deflation 

5 

21.98 

0.333 

8.24 

Inflation 

6.5 

29.26 

0.276 

11.47 

Deflation 

6.5 

29.11 

0.309 

11.12 

Inflation 

8 

34.57 

0.280 

13.50 

Deflation 

8 

34.30 

0.301 

13.18 

Inflation 

9.5 

35.84 

0.298 

13.81 

Deflation 

9.5 

38.98 

0.282 

15.20 

27 


INFLATION 


.  CO 

•  CO 

•  Tf 

•  CM 


»  \V-J 
8  O 

d 


o 


CM 


III 

cc 

3 
CO 
CO 
Ui 

cc 

a.  q 

c 

X 

5  E 

<  o 
X 

o 

3  ' 

5.  o.< 


*  ‘0I1VU  S.NOSSIOd 


28 


Figure  2.2-4.  The  relationship  between  the  Poisson's  ratio  and  the 

tranapulaonary  pressure  (p^  -  pp^)  for  rabbit  lungs  in 
Inflation  and  deflation  processes 


exponential  pseudo -et rain  energy  function  proposed  by  Fung  (1975)  was  used  to 
describe  the  experimental  results*  The  material  constants  obtained  are  com¬ 
pared  with  those  of  the  dog's  lung  given  by  Vavter  et  al.  (1978).  The  vari¬ 
ation  of  these  material  constants  with  age  was  not  significant. 

2.3.1.  Introduction 

For  the  analysis  of  the  stress  and  strain  in  the  lung  subjected  to  static 
and  dynamic  loading,  it  is  necessary  to  know  the  mechanical  properties  of  the 
lung  tissue.  Collection  of  a  basic  set  of  data  on  the  mechanical  properties  is 
a  prerequisite  for  dynamic  analysis.  TWo  types  of  data  exist  in  the  litera¬ 
ture.  In  one.  the  lung  is  first  inflated  uniformly  and  then  a  small  perturba¬ 
tion  is  Imposed  on  it.  In  the  other,  the  lung  is  subjected  to  various  types  of 
finite  deformation.  In  the  former,  the  incremental  stresses  are  linearly 
related  to  the  Incremental  strain;  and  experiments  and  data  reduction  are  sim¬ 
ple.  In  the  latter,  the  stress-strain-history  relationship  is  much  more  com¬ 
plicated.  Depending  on  the  applications,  however,  both  types  of  data  are 
Important . 

In  this  section,  results  from  experiments  on  the  human  lung  are  presented 
in  which  the  lung  tissue  was  subjected  to  large  deformation  in  a  "biaxial" 
loading  condition,  and  data  on  the  stress-strain  relationship,  hysteresis, 
strain  rate  effect,  relaxation,  and  creep  were  obtained.  A  constitutive  equa¬ 
tion  is  presented,  and  the  material  constants  are  identified. 

2.3.2,  Materials  and  Methods 

Human  lungs  were  obtained  at  autopsy  less  than  48  hours  post  mortem  from 
patients  whose  cause  of  death  did  not  produce  an  abnormality  in  the  lung.  Fol¬ 
lowing  excision,  pressure-volume  curves  using  air  were  obtained  for  each  lung. 
Three  of  the  lungs  were  then  degassed  as  follows.  The  main  bronchus  of  each 
lung  was  cannulated  and  connected  to  a  vacuum  source.  The  lung  was  then  hung 
inside  an  airtight  vacuum  chamber  containing  water  and  maintained  at  a  pres¬ 
sure  of  about  -34  cm  ^0  which  corresponds  to  the  saturation  vapor  pressure  of 
water.  The  vacuum  source  connected  to  the  bronchi  was  adjusted  to  a  pressure 
of  1-2  cm  ^0.  Thus  the  transpulmonary  pressure  was  36  cm  H^O.  Each  lung  was 
subsequently  slowly  degassed.  Following  degassing  they  were  submerged  in 
saline  and  allowed  to  be  filled  via  the  bronchus  until  an  equilibrium  configu¬ 
ration  was  reached.  TWo  other  lungs  were  not  degassed,  but  were  allowed  to 
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collapse  and  than  filled  with  aallne  under  a  alight  positive  pressure.  The 
liquid-filled  lungs  were  then  froaen  for  cutting.  Slabs  of  lung  tissue  were 
cut  4-5  an  thick  using  a  bandsaw,  from  which  specimens  about  3*3*  0.4  an  In 
sise  were  obtained. 

Once  cut  for  testing,  each  slab  of  lung  tissue  was  allowed  to  thaw  and 
its  thickness  measured  using  .a  tissue  Micrometer  at  a  minimum  of  five  points 
about  the  specimen  surface.  Only,  those  specimens  of  reasonably  even  thickness 
were  tested.  Specimens  cut  from  lungs  which  had  not  been  degassed  before 
saline  filling  and  freesing  were  subsequently  placed  in  a  vacuum  chamber  for 
1-2  hours  for  degassing  before  testing.  Sight  loops  of  surgical  silk  were  sewn 
to  each  of  the  four  edges  of  the  specimen.  Silk  strings  were  then  used  to  tie 
the  loops  to  the  force  bridges  of  the  experimental  apparatus.  The  apparatus, 
nicknamed  TRIAX,  is  a  servo-controlled  optical-electromechanical  device  de¬ 
scribed  in  detail  by  Vawter.  et  al.  (1978).  In  order  to  avoid  the  edge  effects 
introduced  by  this  hooking  system,  a  rectangular  target  was  painted  in  the 
center  of  the  specimen  using  titanium  dioxide  powder  (with  a  carbon  black 
border  for  contrast)  from  which  displacements  were  monitored  with  video  dimen¬ 
sion  analysers.  This  target  occupied  roughly  7Z  of  the  specimen's  surface 
area.  At  all  times  during  testing  the  specimen  was  immersed  in  a  circulating 
saline  bath  kept  at  37 °C.  The  salinity  was  maintained  by  dripping  distilled 
water  as  needed.  The  pH  of  the  solutions  was  6.3,  7.25  and  7.4. 

Stress-strain  curves  for  slabs  of  human  lung  parenchyma  were  obtained  as 
follows.  Using  the  capability  of  the  TRIAX  to  control  either  the  force  or  the 
displacement  in  x  and  y  directions,  we  maintained  a  constant  force  level  in 
one  direction  and  a  displacement  varying  sinusoidally  with  time  in  the  other 
direction.  The  amplitude  of  the  sine  wave  in  displacement  was  adjusted  to  pro- 
duce  a  maximum  stress  of  60-70  ga/cm  in  the  direction  of  stretch  and  a  mini¬ 
mum  of  about  0.  The  frequency  was  0.4  Hz.  Four  levels  of  constant  force  in  the 
opposite  direction  corresponding  to  5,  10,  20  and  30  gm/cm  stress  were 
chosen,  representing  roughly  7-50Z  of  the  maximum  cyclic  stress  in  the  perpen¬ 
dicular  direction.  These  stresses  are  based  on  the  undeformed  dimensions  of 
the  specimen  and  are  therefore  defined  in  the  Lagrangian  sense.  As  with  virtu¬ 
ally  all  soft  tissues,  lung  parenchyma  must  be  preconditioned  before  each  test 
in  order  to  obtain  reproducible  results.  In  our  testing,  a  steady  state  was 
reached  within  ten  cycles,  with  data  being  collected  on  the  last  cycle  by 
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digital  stapling  using  a  PDP  11/03  computer.  All  loads  on  the  specimen  were 
released  after  each  test  so  as  to  avoid  creep  over  long  periods  of  time.  At 
various  times  the  test  in  which  the  lateral  load  was  held  conrtant  at  a  stress 
level  of  3  gm/cm  was  repeated  so  that  the  repaatlbllity  of  the  measurements 
could  be  assessed.  Anisotropy  is  assessed  by  the  experimental  data,  and  also 
by  repeated  testing  with  the  role  of  x  and  y  Interchanged. 

2.3.3.  Ixpmrlmmmtal  tamulta 

2 

A  typical  stress-strain  curve  for  the  case  ox  “  20  gm/cm  is  shown  in 
Figure  2.3-1.  The  curve  la  highly  nonlinear,  having  shallow  slopes  for  small 
strains  and  much  steeper  slopes  for  large  strains.  A  hysteresis  loop  Is  seen. 
The  stress-strain  relationship  in  loading  la  not  the  same  as  that  in  unload¬ 
ing.  But  we  shall  see  later  (Figure  2.3-5)  that  the  hysteresis  loop  Is  rather 
insensitive  to  the  strain  rate,  and  thus  the  loading  curve  and  unloading  curve 
are  individually  repeatable  In  cyclic  loading.  Henceforth  In  examining  the 
effects  of  other  parameters,  only  the  loading  curves  will  be  considered. 

In  Figure  2.3-2,  the  effect  of  lateral  tension  on  the  stress-strain  rela¬ 
tionship  (loading  curves)  of  a  vs  Xy  is  shown.  The  lateral  load  was  huld  con¬ 
stant  at  5,  10,  20,  30  gm/cm.  It  is  seen  that  the  higher  ox  la,  the  higher 
the  stress  oy  at  any  given  Xy.  The  curves  In  Figure  2.3-3  Illustrate  the  vari¬ 
ation  among  biological  specimens.  Here  the  oy  vs  Xy  curves  for  five  specimens 
when  o  was  fixed  at  20  gm/cm^  are  shown.  Notice  that  X  ranged  from  1.73  to 

*  n  *  2 

1.82  when  oy  was  63  gm/cm  ;  whereas  oy  ranged  over  40  to  62  gm/cm  when  Xy  « 
1,7.  In  Figure  2.3-4  the  lateral  displacements  for  tests  shorn  in  Figure  2.3-2 
are  exhibited.  When  a  larger  lateral  force  is  applied,  the  area  of  the  hyster¬ 
esis  loop  Is  smaller. 

To  assess  the  effects  of  strain  rate  on  the  stress-strain  relationship, 
the  rate  of  stretch  was  varied  over  20-fold  (the  cycle  time  was  varied  from  25 
to  500  s).  The  results  are  shown  In  Figure  2.3-5.  It  Is  seen  that  the  stress- 
strain  relationship  of  human  lung  tissue  is  rather  insensitive  to  strain  rate. 

2.3.4.  Relaxation  and  Creep 

In  Figure  2,3-6  some  typical  relaxation  curves  of  human  lung  tissue  sub¬ 
jected  to  a  step-change  in  strain  are  plotted.  It  is  seen  that  the  relaxation 
phenomenon  is  significant  for  the  lung  tissue,  especially  at  higher  levels  of 
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Figure  2.3-1.  A  typical  etreee-etrain  curve 
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LAGRANGIAN  STRESS  <ry  (gm/cm2) 


STRETCH  RATIO  Xy 


Figure  2.3-2.  Stress-strain  curves  of  oy  vs  X  when  o  was  held  at  various 
fixed  values 
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Figure  2.3-5.  Stress-strain  relationship  oy  vs  Xy  measured  at  various  strain 
rates 
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strain.  Figure  2.3-7  shows  the  results  on  a  logarithmic  time  scale  that  re¬ 
veals  the  various  periods. 

Figure  2.3-8  shows  the  creep  curves  of  the  lung  tissue  under  constant 
loads.  It  is  seen  that  creep  in  the  lung  exists,  but  is  remarkably  small. 

2.3.5.  The  Constitutive  Equation 

A  mathematical  expression  describing  the  mechanical  properties  of  a  mate¬ 
rial  is  called  the  constitutive  equation.  The  mechanical  behavior  of  human 
lung  tissue  exhibited  in  Figures  2.3-1  through  2.3-8  is  fairly  similar  to  that 
of  some  other  soft  tissues  so  that  the  constitutive  equation  suggested  by  Fung 
(1975,  1981)  seems  applicable.  This  constitutive  equation  is  based  on  the 
concept  of  quasl-linear  viscoelasticity  and  pseudo-elastic  strain  energy. 
Expressed  in  the  form  of  a  pseudo  strain  energy  function,  we  have: 


p  W 
o 


7  C  exP[aiEx  +  a2E 


+  2a.  E  E  ] 
A  x  yJ 


symmetric  terms  obtained 
by  permuting  x,  y,  and  z 


(2.3-1) 


where  pQ  is  the  density  of  the  lung  in  the  reference  state,  W  is  the  pseudo 
strain  energy  per  unit  mass,  and  C,  aj,  a2,  a4  are  material  constants.  Ex,  Ey 
are  Green’s  strains  which  are  related  to  the  stretch  ratios  X  and  X  as 

x  y 

follows : 

EX  "  (*!  ■  0 /Z  *  Ey  -  (*y  -  0/2  .  (2.3-2) 

Hence,  ve  have  the  Lagranglan  stresses 

C  X*(alBx  +  a4Ey^  exptalEx  +  a2Ey  +  2a4BxEy 
+  C  Xx(aiEx  +  a4Ez)  exp[aiEx  +  a^2  +  2a4ExKj  ,  (2.3-3a) 

cXy(a2Ey  +  a4Ex)  expfa^2  +  a2E2  +  2a4ExEy 
+  CXy(a2Ey  +  a4Ez)  exP[fllEz  +  a2Ey  +  2aAEzEy]  *  (2-3"3b) 
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Figure  2.3-7.  Relaxation  of  human  lung  tissue 
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Figure  2.3-8.  Creep  of  human  lung  tissue 


STRESS 


Equations  (2.3-3)  show  that  the  constant  C  determines  the  overall  stress 
level,  whereas  a^  and  a2  determine  the  range  of  change  of  stress  with  increas¬ 
ing  stress,  and  a^  determines  the  coupling  between  two  perpendicular  direc¬ 
tions.  If  the  lung  parenchyma  is  assumed  to  be  isotropic  in  the  reference 
state,  then  at  and  a2  are  equal.  Hence,  for  initially  isotropic  lung  tissue, 

«,  -  cM"lEx  +  **V  +  #  +  2*4ExEy] 

+  +  »4Et)  .xpl.jlE2  +  E2)  +  2*4Bxs,]  .  (2.3-4.) 

“y  -  cy*iEy  +  >4E.)  “p[*i(Ex  +  Ey)  +  Wy) 

+  C\  (a.E  +  a.E  )  exp[a.(E2  +  E2)  +  2a.EE]  .  (2.3-4b) 

y '  1  y  4  rL  lv  z  yJ  4  s  yJ 

When  experimental  values  of  o^,  o^,  E^,  E^,  E^  are  substituted  into  Eqs. 
(2.3-3)  and  (2.3-4),  we  can  identify  the  material  constants  C,  ap  a2,  and  a^ 
by  minimising  the  sum  of  the  squares  of  errors  between  experimental  and  com¬ 
puted  stresses. 

2.3.6.  Material  Constanta 

Five  human  lungs  have  been  tested,  as  listed  in  Table  2.3-1.  All  test 
specimens  are  tested  at  37°C.  Two  (1712,  2111)  were  tested  at  pH  »  7.4,  two 

(2508,  1208)  at  pH  -  6.3,  and  one  (2109)  at  pH  -  7.25.  For  each  specimen,  one 

hundred  data  points  were  recorded  in  the  computer.  These  are  used  to  identify 
the  material  constants. 

Our  results  are  presented  in  Tables  2.3-2  and  2.3-3.  The  material  con¬ 
stants  are  listed  separately  for  the  loading  and  unloading  conditions,  and 
under  the  assumptions  of  isotropy  or  anistropy.  In  general,  the  fitting  be¬ 
tween  the  constitutive  equation  and  the  experimental  data  is  very  good. 
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Table  2.3-1.  Identification  of  Hunan  Lung  Spednens 


Table  2.3-2,  Mean  Values  of  Lung  Elasticity  Material  Constants  in 
Loading  Process 


a. 

Lung  Tissue  Considered  Anisotropic 

Specimen 

No. 

T°C 

pH 

C±S.D. 

(g/cm2) 

a^±S.D, 

AjlSeDe 

«4±s.d. 

1712 

37 

7.4 

5 . 84± 1.19 

0.8010.26 

0.7510.12 

0.5310.14 

2111 

37 

7.4 

8.5713.19 

0.6810.24 

0.7110.14 

0.4410.14 

2508 

37 

6.3 

4.57±1.48 

1.60H.54 

1.0010.25 

0.7310.50 

2109 

37 

7.25 

8.7513.46 

0.5410.12 

0.7310.17 

0.3710.07 

1208 

37 

6.3 

10.32±2.90 

0.5310.16 

0.7310.13 

0.4710.11 

Mean: 

7.57±2.33 

0.8310.44 

0.7810.12 

0.5110.14 

b.  Lung  Tissue 

Considered 

Isotropic 

Specimen 

No. 

T°C 

pH 

c±s.p. 

(g/cm2) 

-  a2±S.E 

i. 

a4±S.D. 

1712 

37 

7.4 

5.2110.78 

0.7810.02 

0.5310.07 

2111 

37 

7.4 

5.5210.90 

0.8610.07 

0.5110.05 

2508 

37 

6.3 

3.9212.15 

1.1110.33 

0.6310.17 

2109 

37 

7.25 

3.8014.01 

1.1610.39 

2.6310.19 

1208 

37 

6.3 

3.1312.89 

1.4010.41 

0.9110.24 

Mean: 

4.3211,01 

1.0610.25 

0.6410.16 
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Table  2.3-3.  Mean  Values  of  Lung  Elasticity  Material  Constants  In  Unloading 
Process 


.  Lung  Tissue  Considered  Anisotropic 
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Because  of  the  relatively  slow  speed  of  stress  waves  In  the  lung,  the 
features  of  wave  reflection,  focusing,  and  interaction  often  dominate  the 
scene  of  lung  dynamics  In  response  to  Impact  loading.  In  any  attempt  to 
analyse  the  lung  dynamics,  it  is  necessary  to  know  the  speed  of  wave  propaga¬ 
tion.  In  the  past,  Clemedson  and  JBnsson  (1962)  have  measured  the  wave  speed 
In  the  lung  by  Inserting  two  pressure  transducers  (of  diameter  3/32  Inch  or 
2.38  mm)  Into  the  airways  of  the  right  and  left  lungs  of  an  anesthetlsed 
rabbit  or  calf  and  then  subjected  the  animal  to  an  air  shock  wave.  By  measur¬ 
ing  the  difference  of  the  arrival  times  of  the  stress  waves  at  the  probes  and 
the  distance  between  the  probes  by  autopsy,  they  obtained  the  velocity  of  the 
main  part  of  the  pressure  pulse  In  the  calf  lungs  inflated  to  approximately 
lntravltal  size  to  be  between  24  and  30  m/sec.  For  rabbit  lungs  they  found 
speeds  of  15,  32,  and  69  m/sec  In  several  animals  but  considered  the  data  as 
uncertain.  In  collapsed  rabbit  lung  they  found  a  wave  speed  of  15  m/sec.  Rice 
(1983),  using  a  microphone  and  a  spark  gap  sound  source,  measured  the  transit 
time  of  sound  through  the  lung  and  estimated  the  sound  speed  In  air-filled 
excised  horse  lung  to  be  between  25  and  70  m/sec,  depending  on  lung  volume. 
Dunn  and  Fry  (1961),  using  reflection  coefficient  for  ultrasound  energy  and 
the  density  of  the  lung  tissue,  calculated  a  wave  speed  of  650  m/sec.  These 
are  the  data  In  the  literature  as  far  as  we  know. 

The  lung  has  a  complex  structure  and  the  lung  parenchyma  is  a  composite 
material  with  a  highly  nonlinear  stress-strain  relationship.  Many  types  of 
waves  can  exist  in  it.  Which  one  is  excited  depends  on  the  mode  of  excitation 
and  the  amplitude  of  the  loading.  The  10-  to  20-fold  difference  between  the 
data  quoted  above  reflects  the  fact  that  different  kinds  of  waves  were  excited 
in  different  experiments.  It  follows  that  one  should  design  one's  experiment 
according  to  the  intended  application,  so  that  the  type  of  the  impact  load  and 
its  amplitude  level  can  be  properly  chosen.  In  the  present  study,  we  are 
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interested  in  the  dynamics  of  the  lung  in  response  to  iapsct  or  shock  wavs 
acting  on  the  chest;  hence  the  iapsct  load  level  should  be  relevent  to  the 
problea  of  possible  lung  Injury.  In  the  following  we  present  the  results  of 
our  aeasureaent  of  wave  speed  in  the  lungs  of  the  rabbit*  the  goat*  and  aan 
with  a  aethod  in  which  a  water  jet  was  used  to  lapose  an  iapact  load  on  the 
lung.  A  pressure  transducer  was  used  to  record  the  pressure  wave  produced  by 
the  jet  hitting  the  lung.  The  transducer  was  also  placed  at  a  point  on  the  far 
side  of  the  lung  to  sense  the  arrival  of  the  stress  waves.  A  photosensor  was 
placed  at  various  places  between  the  nossle  of  the  water  jet  and  the  lung  to 
aeasure  the  arrival  of  the  water.  Froa  these  signals  the  velocity  of  the  water 
jet  and  the  speed  of  the  stress  wave  in  the  lung  can  be  deterained. 

The  objective  of  our  study  is  siailar  to  that  of  Cleaedson  and  J&nsson, 
but  different  from  those  of  Dunn  and  Fry  and  those  of  Rice.  The  impact  pres-* 
sure  level  on  the  surface  should  be  in  the  range  of  1-5  psi  (7-35  kPa),  and 
the  surface  velocity  should  be  in  the  order  of  1-10  a/ sec.  According  to 
Clemedson  and  JSnsson  (1975)  severe  hemorrhages  in  the  lung  are  highly  proba¬ 
ble  for  a  maximum  surface  velocity  exceeding  about  10  a/ sec.  We  chose  to  use  a 
water  jet  to  deliver  an  Impact  in  this  range.  The  velocity*  aoaentua  and  area 
of  impact  of  the  water  jet  can  be  controlled. 

Our  choice  of  water  jet  is  based  on  the  following  consideration.  Cleaed¬ 
son  and  Jbnsson  (1964)  used  the  detonation  of  a  high  explosive  in  a  "wave 
guide"  (constituted  of  two  concentric  cylinders)  to  lapose  a  planar  shock  wave 
on  the  animal.  By  their  aethod  the  determination  of  the  wave  arrival  time  is 
difficult  as  can  be  seen  from  the  cathode  ray  traces  in  their  paper.  In  earli¬ 
er  stages  of  our  own  experiment*  we  designed  and  used  several  shock  tubes 
(using  compressed  air  and  paper  membrane)  to  produce  a  loading  in  an  over¬ 
pressure  range  of  a  few  psi.  The  weak  supersonic  shock  is  followed  by  a  bulk 

w 

flow.  The  lung  responds  to  both  the  weak  supersonic  shock  wave  end  the  large 
bulk  flow,  and  the  determination  of  wave  propagation  characteristics  is  often 
difficult.  This  difficulty  was  minimized  by  using  the  water  jet. 


3.2  APPARATUS 


Figure  3-1  is  a  photograph  of  the  test  equipment  used  for  measuring  the 
elastic  wave  propagation  in  the  lung.  Figure  3-2  is  a  block  diagram  of  the 
testing  apparatus.  The  apparatus  consists  of  three  subgroups,  the  first  of 
which  is  the  water  jet  for  imposing  impact  loading  on  the  surface  of  the  lung. 
To  obtain  a  well  defined  Impulse  loading  covering  a  circle  of  2  mm  diameter  on 
the  surface  of  the  lung,  a  water  jet  was  constructed.  A  fluid  reservoir  driven 
by  a  pressure  regulated  air  source  charged  a  small  hydraulic  accumulator 
through  a  normally  open  flexure  electronic  poppet  valve.  This  valve,  manufac¬ 
tured  by  Clippard  Minimat ic,  has  an  activating  time  of  <  3  ms  due  to  a  very 
low  mass  of  flexure  poppet.  The  valve  was  designed  to  operate  at  5  volte  DC  at 
100  mA.  When  a  normal  trigger  was  applied  to  the  valve  sequencer,  valve  A 
would  close,  isolating  the  accumulator  from  the  fluid  reservoir.  After  30  ms 
the  valve  B,  an  identical  flexure  poppet  valve,  would  open  for  5  ms,  discharg¬ 
ing  the  accumulator  through  the  nozzle.  The  nozzle  itself  consisted  of  5  holes 
of  0.3  mn  diameter  each  grouped  in  a  2  m  circle. 

The  second  subgroup  consisted  of  the  photosensor  system  which  was  used  to 
sense  the  arrival  of  the  water  jet.  The  photosensor  could  be  moved  with  preci¬ 
sion  along  the  axis  of  the  jet  to  calculate  the  velocity  of  the  fluid  in  the 
air.  The  conditioning  electronics  incorporated  a  hi -pass  filter  with  a  time 
constant  of  100  ys  to  eliminate  ambient  light  effects. 

The  third  subgroup  consisted  of  the  pressure  transducer  used  to  record 
the  response  of  the  lung  due  to  water  impact.  The  arrival  time  and  the  pres¬ 
sure  wave  were  recorded.  The  pressure  transducer  was  connected  to  a  Neff 
amplifier  and  the  output  was  shown  on  a  Tektronix  5440  oscilloscope.  The 
pressure  transducer  was  an  Endevco  Model  8510-2  m21  with  a  calibrated  sensi¬ 
tivity  of  151.6  m  volts/psi. 

3.3  ANIMAL  PREPARATION 

Animals  such  as  goats,  etc.  were  first  anesthetized  using  a  proper  dose 
(~  7  mg/kg)  of  Nembutal.  Heparin  (145  units/kg)  was  then  injected  into  the 
bloodetream  to  prevent  blood  coagulation.  The  trachea  was  cannulated.  Through 
a  midline  incision  of  the  chest,  the  lung  was  exposed  and  excised.  An  isola¬ 
tion  of  the  lung  lobes  was  then  made  after  placing  the  lung  in  a  saline  filled 
trey. 
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Figure  3-1.  The  test  equipment  for  measuring  the  elastic  wave  propagation 
In  the  lung 


ALVE 


The  experimental  protocol  is  shown  in  Figure  3-2.  The  lung  was  placed 
firmly  on  a  pressure  transducer  in  a  supporting  basin.  A  water  jet  was  placed 
directly  above  the  lung,  with  a  distance  of  1.8  cm  from  the  bottom  of  the 


nozzle  to  the  surface  of  the  lung.  The  nozzle  has  five  holes  to  allow  five 


small  jets  converging  onto  one  area  on  the  lung  surface.  A  photosensing  tube 
was  placed  between  the  lung  surface  and  the  nozzle.  In  each  experiment,  the 
reservoir  air  pressure  was  kept  constant  (50-60  psi).  With  a  touch  of  the 
control  unit  key,  a  water  jet  was  fired.  The  screen  of  the  oscilloscope  showed 
the  output.  For  each  constant  reservoir  air  pressure,  at  least  three  photo¬ 
sensor  positions  were  used  in  the  experiment.  Then  another  reservoir  air 
pressure  was  selected  and  the  experiments  were  repeated.  The  airway  pressures 
of  the  lung  were  chosen  to  be  0,  5,  10,  15  and  20  cm  l^O.  The  pleural  pressure 
was  zero  (atmospheric).  We  then  plotted  on  an  x-y  coordinate  paper  (Figure 
^  3-3)  the  distance  from  the  top  of  the  pressure  transducer  (d)  versus  the  dif- 


If  France  in  time  (t)  between  the  leading  edges  of  the  pressure  wave  and  photo¬ 


sensor  wave.  For  each  water  ejection  setting  (pressure  in  the  accumulator)  a 
straight  line  was  obtained.  By  changing  the  setting,  different  straight  lines 
were  obtained.  These  lines  would  intersect  at  one  point  which  marked  the  time 
and  distance  of  the  point  where  the  water  jet  hit  the  lung  surface.  In  this 
way  we  obtained  the  lung  thickness  d.  The  corresponding  time  t  is  the  travel¬ 
ing  time  of  impact  wave  in  the  lung  tissue.  Hence  the  wave  speed  of  the  lung 
tissue  is  c  -  d/t. 


3.5  DATA  AHALTSIS 


Figure  3-4  shows  a  typical  set  of  experimental  curves.  Figure  3-4(a) 
shows  the  pressure  wave  due  to  impact  of  the  water  jet  on  the  transducer  lo¬ 
cated  at  the  place  where  the  jet  hits  the  lung  (with  the  lung  removed).  Thi9 
particular  wave  shows  a  peak  pressure  of  1.7  psi.  Figure  3-4 (b)  yhows  the 
usual  format  of  the  outputs  recorded  on  the  screen  of  the  oscilloscope.  The 
upper  curve  is  the  response  of  the  photosensor  indicating  the  arrival  of  the 
water  jet.  The  lower  curve  is  the  response  of  the  pressure  transducer  which 
indicates  the  stress  wave  in  the  lung  at  the  place  whore  the  transducer  was 
located.  The  first  sharp  wave  is  used  to  correlate  '  impact  wave  imposed  on 
the  lung.  The  time  difference  between  the  arrival  of  the  water  jet  on  the 
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(a)  The  pressure  transducer  response  due  to  Impact  of  the  water  jet. 
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(b)  Example  of  output  response  recorded  on  the  screen  of  the  oscilloscope. 
The  upper  curve  Indicates  the  arrival  of  water  jet  on  the  photosensor. 
The  lower  curve  is  the  response  of  the  pressure  transducer  which 


Indicates  the  stress  wave  propagation.  The  first  sharp  wave  is  used 
to  correlate  the  Impact  wave  imposed  on  the  lung.  In  this  particular 
example,  the  peak  is  0.003  psl  according  to  the  calibrated  scale 
(one  division  is  0.0015  psi). 

Figure  3-4.  Typical  experimental  curves 

52 


photosensor  and  first  arrival  of  the  wave  by  impact  was  measured  and  desig¬ 
nated  as  tj  corresponding  to  photosensor  position  d^;  at  least  two  other 
photosensor  positions  were  chosen  and  the  results  are  designated  as  t2»  ^  and 
tg>  dg.  These  data  are  used  to  obtain  the  water  jet  speed  in  the  air  (the 
slope  of  the  curve)  as  well  as  the  distance  and  traveling  time  in  the  lung 
tissue  (Fig.  3-3). 

3.6  RESULTS 

Table  3-1  gives  the  experimental  results  of  goat  lung  for  several  airway 
pressures.  The  values  of  wave  speed  c  lie  in  the  range  of  30  to  70  m/s,  and 
have  a  general  pattern  of  increasing  with  the  airway  pressure.  Here  n  is  the 
number  of  experiments  and  a  is  the  standard  deviation. 

Table  3-2  summarizes  the  experimental  results  of  the  rabbit  lung. 

3.7  COMPARISON  WITH  DATA  IN  THE  LITERATURE 

It  is  often  stated  that  the  velocity  of  stress  wave  in  soft  tissues  of 
the  body  is  equal  to  the  velocity  of  sound  in  water  due  to  the  high  water  con¬ 
tent  of  most  animal  tissues  [see,  e.g.,  Corey  (1946),  Schardin  (1950);  c.f. 
also  White  and  Richmond  (1959)].  The  data  on  the  lung  show  that  this  is  in¬ 
correct  . 

Table  3-3  gives  the  velocity  of  sound  in  various  soft  and  hard  tissues. 
Comparing  the  velocity  of  sound  in  the  lung  with  that  in  other  tissues,  one 
sees  at  once  that  the  wave  speed  is  uniquely  low  in  the  lung. 


3.8  DISCUSSION 

The  wave  speed  obtained  from  lungs  with  their  blood  drained  were  differ¬ 
ent  from  those  obtained  for  lungs  filled  with  blood.  Usually  the  values  for 
the  blood-drained  lungs  were  found  to  be  higher.  In  one  experiment  using  goat 
lung  with  no  blood  in  the  pulmonary  blood  vessels,  the  obtained  values  corre¬ 
sponding  to  airway  pressures  of  0,  5,  10,  and  15  cm  ^0  were  38.7  ±1  .8,  52.5 
±3.0,  67.7  ±  4.9,  and  77.8  ±3.7  m/sec,  respectively.  Compared  with  the  previ¬ 
ous  results  with  blood  in  the  puli*onary  blood  vessels,  it  is  seen  that  these 
values  are  higher.  This  finding  appears  to  be  reasonable.  If  we  consider  the 
lung  tissue  to  be  incrementally  linearly  elastic,  then  the  square  of  the  wave 
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Table  3-1 .  Stress  Wave  Speed  In  Goat  Lung  Is  Airway  Pressure  and 
ppL  is  Pleural  Pressure  Which  is  Atmospheric 


Mean  Value 

Transpulmonary 

of  Wave 

Standard 

Pressure 

Velocity 

Velocity 

Deviation 

Pa  "  ppl 

(m/s) 

n 

(m/s) 

0 

31.6 

0  30.8  3  31.4  0.4 

31.8 


5 


10 


15 


32.7 

33.8 
31.4 
37.6 


33.9 


2.3 


37.7 
37.4 

36.1 

32.1 

37.8 
35.6 


6  36.1 


1.9 


46.9 

43.0 

46.7 

48.7 

47.2 

48.2 


6  46.8 


1.8 


66.5 
64.7 

58.5 
69.1 


I 


20 


4 


64.7 


3.9 


Table  3-2.  Stress  Wave  Speed  In  Rabbit  Lung  is  Airway  Pressure  and 
ppL  is  Pleural  Pressure  Which  is  Atmospheric 


Mean  Value 

Transpulmonary 

of  Wave 

Standard 

Pressure 

Velocity 

Velocity 

Deviation 

PA  "  PPL 

(m/s) 

n 

(m/s) 

0 

15.2 

15.4 
21.7 
15.0 
15.0 

16.4 


16.5 


2.4 


32.1 

31.1 
28.6 
30.9 
28.6 

22.1 


28.9 


3.3 


12 


16 


29.6 
31.9 

31.3 

31.6 

32.3 

35.5 

35.5 

36.2 

34.2 

39.5 

35.2 

38.6 

34.2 
37.0 

37.2 
36.8 


31.3 


0.9 


35.3 


0.8 


36.9 


1.7 
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Table  3-3.  Velocity  of  Sound  in  Varioue  Tissues,  Air  and  Water 


Tissue 

Velocity 
of  Sound 
(m/sec) 

Density 

(g/cm3) 

Reference 

Muscle 

1580 

1 

Ludwig  (1950),  Frucht 
(1953),  von  Gierke  (1964) 

Fat 

1450 

1 

Ludwig  (1950),  Fruch  (1953) 

Bone 

3500 

2.0 

Clemedson  &  JBnsson  (1961) 

Ribs  &  intercostal 
muscle 

<1000 

Clemedson  &  Jbnsson  (1961) 

Collapsed  lung 

650 

(ultrasound) 

0.4 

Dunn  &  Fry  (1961) 

Collapsed  lung 
pneumonitis 

320 

(ultrasound) 

0.8 

Dunn  &  Fry  (1961) 

Lung,  air  filled,  horse 

25 

0.6 

Rice  (1983) 

Lung,  air  filled,  horse 

70 

0.125 

Rice  (1983) 

Lung,  air  filled,  calf 

24-30 

Clemedson  &  Jttnsson  (1962) 

Air 

340 

Dunn  &  Fry  (1961) 

Water,  distilled,  0°C 

1407 

Kaye  &  Laby  (1960) 

Air  bubbles  (45%  by 
vol.)  In  glycerol  &  H20 

20 

Campbell  &  Pitcher  (1958) 

r 
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2 

speed  c  is  proportional  to  the  Young's  modulus  E  and  inversely  proportional 
to  tissue  density  p.  Both  E  and  p  change  with  the  tranamural  pressure,  but 
when  the  lung  is  drained  of  blood  the  primary  effect  is  a  reduction  of  the 
density  of  the  lung  parenchyma;  and  a  rise  in  wave  speed  is  expected.  Rice 
<1983)  used  He,  air,  and  SFg  in  his  experiments,  and  showed  that  the  lower  the 
density  of  the  gas  the  higher  the  wave  speed. 
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